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A new mechanism is proposed for the magnetization reversal of molecular nanomagnets such as Peg. In 
this process the spin tunnels from the lowest state near one easy direction to the first excited state near the 
opposite easy direction, and subsequently decays to the second easy direction with the emission of a phonon, 
or it first emits a phonon and then tunnels to the final state. This mechanism is the simplest imaginable one 
that allows magnetization relaxation in the presence of a longitudinal magnetic field that is so large that the 
nuclear spin environment cannot absorb the energy required for energy conservation to hold. It is proposed 
as a way of understanding both magnetization realaxation and Landau-Zener-Stiickelberg experiments. The 
requisite Fermi golden rule rate, and the spin-flip rates are calculated, and it is found that these rates are much 
too low by several orders of magnitude. Thus the understanding of magnetic relaxation in the experiments 
remains an open question. 
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1 Introduction and Background 

Over the last fifteen years or so, molecular nanomagnets (also known as single-molecule magnets, or 
molecular magnets) have provided us with an entirely new class of magnetic systems, showing many 
novel phenomena not seen previously Among the most dramatic of these is the observation of gap 
oscillations wherein the tunnel splitting between the two lowest energy states on opposite sides of an energy 
barrier oscillates as a function of a static magnetic field applied along the hard axis of the molecule ||2l- 
Although many aspects of this phenomenon can be understood by considering the idealized problem of 
an isolated molecule, a full understanding of the experimental procedures and results presents several 
complexities and puzzles. It is the purpose of this paper to address one of these puzzles. To keep the 
discussion focussed and free of lapidary generalizations, we will consider the example of the Fes molecule 
in this paper, but the central ideas are applicable more broadly. 

To understand the basic phenomenon at issue, let us first consider one Feg molecule in the solid, and 
ignore its interaction with other molecules and any other environmental degree of freedom. The molecule 
has spin 10 in its ground manifold. At low temperatures only the spin degree of freedom has any life in it, 
and all others are frozen. This degree of freedom is governed by an anisotropy Hamiltonian 

Hs = kiS'^ + k2Sy + ^4 - .g/isS • H, (1) 

where S is a dimensionless spin operator (of magnitude S, equal to 10 for Fes), ki > k2 > are anisotropy 
coefficients, g is a g-factor (equal to 2 for Fes), and H is an external magnetic field. The term H4 is of 
fourth order in the components of S, and while it is responsible for surprising and important features in the 
tunneling spectrum, we do not show it explicitly because we will not need to discuss those aspects in this 
paper. 

When H = 0, Eq. ([T]) has two ground states corresponding (classically speaking) to Sj| ±z, or m = ±10, 
where m is the Sz or Zeeman quantum number These two states are degenerate, and mixed by tunneling, 
as shown in Fig. [T] If we now turn on H||x, the new classical ground states are still degenerate, but 
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Fig. 1 Scliematic energy level diagram of Feg showing (a) the basic tunneling process between ground 
levels, and (b) the phonoemissive process, in which the molecule tunnels from m = —10 to an intermediate 
virtual state m = 9, and then makes a transition to m = 10 with the emission of a phonon (as shown 
by the dot-dashed line). Or, the phonon can be emiited first accompanying a transition to the m = —9 
state, followed by tunneling to the m = 10 state. This process appears to be the simplest way in which the 
magnetization can relax when the bias e is much greater than W , the effective width of the levels induced 
by the nuclear spin environment. Direct tunneling in this case would not conserve energy, and the nuclear 
spins are incapable of absorbing an energy much greater than W . 



the tunnel splitting between them (A_io,io = A) does not increase monotonically with H^. Instead it 
oscillates as shown in Fig.|2] This oscillation is best understood in terms of instanton ||3]|4]|5l, but readers 
unfamiliar with instantons may regard it as an outcome of direct diagnolization of the Hamiltonian ([Til. 
Similarly, if a magnetic field is applied along the z axis so as to bring the m = — 10 and 7ti = 9 states 
into resonance, the system will be able to tunnel between these states. This tunnel splitting (A-io g = A') 
also oscillates as a function of an additional field along the hard axis. The same holds for tunneling 
between m = —10 and m = 8, and other pairs of states. By symmetry, if the z axis field is such to bring 
the m = —9 and m, = 10 states into resonance, the spUtting A_9.io will be the same as A-io.g. 
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Fig. 2 The gap oscillations in Feg. Part (a), kindly supplied by Dr. Wemsdorfer, shows data from Ref. (2). The curve 
labeled ti = is the — 10 10 tunneling, that labeled n = 1 is the — 10 — 9 tunneling, and so on. Part (b) shows a 
fit to data for the tunnel splitting between the lowest pair of levels, taken from Ref. 1121 . 



(To avoid misunderstanding, we note that when we speak of the state with m = —10, say, we do not 
literally mean the eigenstate of Sz with eigenvalue —10. Rather, we mean the eigenstate of Hs (ignoring 
tunneling) that maps on to the Sz = —10 state in the following way. We imagine subtracting a zeroth order 
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Hamiltoian (fci + k2){S'^ + Sy)/2 from Eq. ([T]), and treating the result as a perturbation. The state that we 
are talking about is then the one that would develop from the Sz — —10 state by low order perturbation 
theory. It would perhaps be better to use a notation such as m* for the perturbed state, but as long as this 
qualification is understood, there is little benefit from doing so.) 

If an Peg molecule really did not interact with the rest of the world, the tunneling described above would 
lead to coherent flip-flop as in the inversion resonance of ammonia. No such flip-flop is seen, and indeed 
that is to be expected. Any environmental degre of freedom which couples to the magnetic moment of 
the molecule will tend to suppress quantum coherence, and for resonance between states with such a large 
difference in their magnetic moments, one would expect that all vestiges of coherence are destroyed. This 
is indeed so, and one finds that because of the nuclear spin environment, transitions between the ±5* states 
are totally incoherent, and one finds a transition probability per unit time given by ||6l[3 



4 W \2W^ 

Here, W ~ lOEdn, where E^n is the energy of dipole-dipole interaction between the molecular electronic 
spin and the nuclear spins of nearby nonmagnetic atoms such as N and H which are always present in the 
molecules studied, and e is the bias, or the energy of the m ^ —S state relative to the m = S state. This 
bias can arise from an externally applied field along the z axis, or from the dipole field created by other 
molecules in the sample. Indeed, it is found that this dipole field is of order 100 Oe |[8l|9l, so the bias is of 
order 0.1 K. By comparison, Edn ~ 1 mK. 

Given that e ^ W for most molecules in any solid sample, the spins of most molecules are frozen, 
in that they do not even undergo incoherent tunneling. To overcome this problem, the authors of Ref. 121 
use the ingenious idea of sweeping through the resonance by applying a time-dependent longitudinal (z 
axis) magnetic field as shown in Fig. [3] This is the Landau-Zener-Stiickelberg (LZS) protocol fTOl . The 
molecule flips from m = —10 to m = 10 (or vice versa) with a probability 

/ 7rA2\ 

PLZS = l-exp(^-^j (3) 

for every sweep through the crossing. Here e is the rate at which the bias changes. For transitions between 
m = ±5, |e| = 4:ijLBS\dHz/dt\, since g = 2. It is found that for between ~ 3mT/s and 1 T/s, the LZS 
formula is obeyed, allowing one to extract A from a measurment of plzs- 

OT = 10 \ ^ m = 10 

+ phonon 




m = —10 m = —10 




(a) (b) 

Fig. 3 The Landau-Zener-Stiickelberg protocol. In part (a) we show the direct LZS mechanism, in which 
the energy levels of the m = ±10 states are swept by applying a time-dependent field along the z axis. In 
(b) we show the same process when the final state also contains a phonon. 

At this point it should be stated that the LZS formula ^ only holds for coherent transitions. For 
incoherent transitions, a different formula, due to Kayanuma ifTTI is better: 
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Remarkably, for fast sweeps, i.e., when |e| ^ jh, both formulas agree: 



7rA2 



Plzs = = 



2?i|er 



(5) 



Experimentally, the LZS formula starts to fail for dH. /dt < i mT/s (the Kayanuma formula is not pre- 
cisely obeyed either), but if one only employs it for fast sweeps, then the details of the decoherence do not 
matter, and the extraction of A is reliable. 

2 The Phonoemssive Tunneling Process 

We are now ready to describe the puzzle mentioned at the start of Sec. [T] Theoretically, at the minima, 
A vanishes strictly. Experimentally, this is not so. In Ref. |12i, this was explained by arguing that every 
molecule sat in a nonzero transverse magnetic field created by the other molecules that acts in addition 
to the applied H^- Thus A 7^ for any molecule, and the measured A should really be found by taking 
an average over the inhomogeneous spread in of order 100 Oe. (Actually, the quantity that should be 
averaged is A^.) In addition, the experimental data show that the inferred values of A at the minima grow 
linearly with H^. To explain this, it was assumed that there was a small degree of misalignment of the 
magnetic axes of the molecules (of order ag ^ 1-2°) because of various types of defects in the material, 
which is rather soft and organic. This means that there is a small nonzero field iJ,ncd — H^crg along the 
local medium axis of each molecule, and since for an isolated molecule A cx iJnicd in the vicinity of one 
of the minima, the measured A would pick up a contribution linear in Hx- 

As can be seen from Fig. |2tb), the fit between the theory based on these assumptions and the experi- 
mental data is good; indeed it is rather too good. For while the first assumption — of an inhomogeneous 
transverse dipolar field — is well justified and hardly an assumption, the second — of misalignment of the 
magnetic axes — is somewhat ad hoc. We have therefore sought another way in which the spin could flip 
for which the rate has an dependence. The mechanism we propose is shown in Fig. [TJb). It involves 
a virtual tunneling transition from the m = —10 to the m = 9 state, followed by a transition from the 
m = 9 state to the m = 10 state with the emission of a phonon. Alternatively, we could first make a virtual 
transition to the m = 9 state accompanied by the emission of a phonon, and then tunnel to the to = 10 
state. Either way, the final state is the same, i.e., the spin is in the m = 10 state, and a phonon has been 
emitted. We refer to this as phonoemissive tunneling. The reverse process, involving the absorption of a 
phonon, will also take place if fegT > e, for then phonons of energy e will be available to be absorbed. In 
this paper we work only at T = 0, so only the emission process takes place. The generalization to T 7^ 
is straightforward. 

At this point it may be useful to clarify that the —10^9 (or —9 — )• 10) tunneling transition is possible 
even if the levels are not in resonance. This transition would have very small probability if we were thinking 
of the TO = 9 state as a true final state (not to mention that it would require some environmental degree of 
freedomn to supply the energy necessary for energy conservation), but we are only exploiting to = 9 as a 
virtual intermediate state. For the virtual process, there is always a matrix element of the Hamiltonian, or 
tunneling amplitude, equal to — iA'/2?i per unit time. Indeed, for Fcg, A' as a function of is maximal 
when A is minimal. 

The process just described could also be germane to experimental situations other than LZS. Consider, 
for example, magnetic relaxation without any swept field. If, for a given molecule, e ^ for some reason, 
the rate Q is very small, and phonoemisive tunneling may be important. Indeed, we shall calculate the 
second order Fermi golden rule rate for this situation first, in Sec. [3] before considering the LZS protocol 
in Sec. m The process could also be important for ultraslow LZS sweeps, and for inverse LZS sweeps lfT3l . 

The second step of phonoemissive tunneling is governed by the spin-phonon Hamiltonian, 
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Here, u°j is the operator for the strain tensor, Dijki is the magnetoelastic tensor, {, } denotes the anticom- 
mutator, and there is an implicit sum over the Cartesian indices i, j, k, and I. The tensor Dijki is not well 
known for Feg, so we shall simplify Eq. (|7]i to 

Hsp ^Du°f{S,,Sj}. (7) 

The quantity D has dimensions of energy, and may be taken to be of the same order of magnitude as the 
anisotropy coefficients fci and fc2- The strain tensor is given by 

where u°p(x) is the displacement field operator given by 



a 

Here, Ua and a^^ are the destruction and creation operators for phonons of mode a, qc Wq, and ec, are the 
wavevector, frequency, and polarization vector for this mode, M is the mass of all atoms in a unit cell, and 
N is the number of unit cells in the crystal. In writing Eq. (|9|, we have kept only long wavelength phonon 
modes, for which relative motion of the atoms in a unit cell is negligible, and this is why the total unit cell 
mass M appears. These are the only modes relevant for low temperatures and low energy processes. We 
can also take the polarization vectors to be real for them. 

For the m = 9 to m = 10 transition, the only nonzero matrix elements arise from the spin operators 
{S^, 5J and {Sy,S^}. We have 

{m = S\{S^,S,}\m = S-l) = {S - ^)V2S = A, 

{m = S\{Sy,S,}\m = S -I) = -i{S - ^)V2S = -iA. (10) 

It is now apparent that we could also take an m = 8 intermediate state, since the operator {Sx, Sy} would 
yield a non zero matrix element, and there is a nonzero amplitude to tunnel from m = —10 to 771 = 8. 
We will not give formulas for this case explicitly, as the requisite modifications are straightforward. Other 
intermediate states are much less important, since they would involve multiphonon processes with far 
smaller probabilities. 

3 Fermi Golden Rule 

We now calculate the rate for phonoemissive tunnelling when there is no swept field. We denote the initial, 
intermediate (or virtual), and final states by \i), \v), and la), where, 

\i) ~ |77i = —10, no phonons), (11) 
\v) = 1 777 = 9, no phonons), (12) 
|a) = |777 = 10, one phonon in mode a). (13) 

The energies of these states are 

Ei ~ £-10, Ey = £9, Ea = £10 + ^l^a, (14) 

where £„ is the energy of the molecular spin alone. Alternatively, the transition could take place through 
a different virtual intermediate state. 



\v) = 1 777 = —9, one phonon in mode a), 



(15) 
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in which case, 

Ev = £-9 + huJa- (16) 

The transition rate is given by the second order Fermi golden rule, 

r^ot = YJ2\v!.?\'m-E^), (17) 

where 

T.(2) {a\nsp\v){v\ns\i) 

K^ • (18) 

It should be noted that the sum over virtual states in Eq. ( fTSl l does not entail a sum over the phonon modes; 
rather it is performed for a particular mode a. The sum over phonon modes is performed in Eq. ( [TtI i. 

It is evident that because of symmetry, the two contributions to Eq. ( fTSl l differ only in the energy 
denominators. We shall assume that the bias e is small in comparison with the energy difference 

E,^ = eg - eio. (19) 

In that case, both energy denominators may be replaced by — i?cx- We need only do the calculation for 
the intermediate state ( fT2] i. and double the answer to get V^'^h Further, it is more accurate to calculate the 
energy difference using the results for e,„ when = 0. 

The actual calculation is straightforward. Consider the intermediate state ( fT2b . For the matrix element 
of Hs we have 

{v\ns\i) = A' /2. (20) 
Next let us examine the matrix element of Hsp ■ Consider the strain field u°p . We have 



<^ = 2 ^ V 2MNio ~ ai)iea,xqa.z + ea.zqa,x)- (21) 
When we take the matrix element of this operator, we are left with the c-number 



Similarly, from it°P we get Uy^, which is the same expresssion with the index a; replaced by y everywhere. 
The matrix elements of the spin parts of Hsp have already been found in Eq. dTol i. Hence, doubling the 
result as explained above to account for the two intermediate states, we get 

l^if =-^«.-<J- (23) 



Therefore, 



rpot-^f^^) ^|«,-*<JP<5(?ic.„-i?./), (24) 



where 



Eif — e_io — eio — £• 



(25) 
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Since Wq, > 0, the rate is nonzero only if e > 0, i.e., if energy conservation requires the emission of 
phonons. To perform the sum over phonon modes, we replace the mode index a by the pair (q, s), where 
q is the wavevector, and s (= 1, 2, or 3) labels the three accoustic modes. The sum over q can be turned 
into an integral in the usual way, and we get 

2tt f ADA'Y hvo ^ f d?q I . n2ic. 

(26) 

where vq is the volume of a unit cell. It should be remembered that Gs depends on q. 

To proceed further and obtain an order of magnitude estimate, we make the simplifying assumption that 
the material is isotropic, so that for any q, there is one longitudinal mode and two (degenerate) transverse 
modes, and that cjqs is either c^g or cxq, where cl and ct are the longitudinal and transverse sound 
velocities. The integral and sum over q and s are then elementary, and we get 



Tpct - TZ — + ^ -74' (27) 



where 

P = M/v^ (28) 

is the mass density of the material. The depedence of the rate ( l27T i is chacteristic of other rates involving 
phonon emission, and the ^ factor has the same origin as in Ref lfT4l . Since ct ^ cl/2, the 
transverse sound term is likely to dominate in Eq. dZTl i. 



4 The Landau-Zener-StUckelberg Rate 



Next, let us consider how the phonoemissive process affects the spin-flip probability if the longitudinal 
field is swept in a LZS protocol. The role of A in Eq. (|2]i is now played by the second-order matrix 
element V^^ . We will only consider the limit of fast sweep, in which case the answer may be found by 

(2) 

a perturbative expansion in V^^ . At the same time, it should be noted that the sweep is slow on the time 
scale h/Ecx required to establish the underlying second-order process. 

With the same notation as in Sec. [3] the time-dependent Hamiltonian for the LZS protocol may be 
written as 



(29) 



where 



Hence, 



Ei(t) = e-w + g^isSH^t, 
Ea{t) ^ eio + huja - gp-sSHzt- 



(30) 
(31) 



(32) 



Further, the energies e±io of a single molecule are now the eigenvalues of Hs with H = 0. It follows that 
that these energies are equal, and we may take them as our reference level: 



eio = £-10 = 0. 



(33) 
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In an interaction picture, the time-dependent state of the system can be written as 

a 

Feeding this form into the Schrodinger equation, we obtain the equations of motion for the ampUtudes a{t) 
and ba{t): 



-iliUat-et^ /2h] 



6„ = _lv^j2)^(i)e*K*-^tV2?«]. (35) 

We wish to solve these equations with the initial conditions 

|a(-oo)| = l, 6a(-oo) = 0. (36) 

Following Kayanuma ifTTI . the solution may be obtained as a power series in V^^'^ by iteratively substituting 
Eqns. ( |35] ) into one another. For us it suffices to take only the first order answer, so we put a{t) = 1 in the 
second equation, which can then be integrated to yield 

h^{t)^^V^f f dfie^["°*i-*?/2''l. (37) 



Thus, ignoring a unimodular multiplicative factor, 



^.M = 4eV^. (38) 
The net spin-flip probability is 

Plzs==EI^"MP (39) 

The sum over phonon modes is performed in the same way as in Sec. |3] except that since there is no delta 
function of energy in the summand, one must cut the sum off in some way. The natural cutoff is provided 
by the peak value of the bias between the m = ±10 states due to the swing in the longitudinal field. 
Writing this swing as i/ac, we define 

Sac = gfisSH,,. (40) 

We then find 

The indirect flip probability must be added to the direct process result (|5]). The result can then be cast 
in the same form as Eq. (|5]l itself, provided we replace A by an effective splitting Aoff , where 



AL = A2 + A'2 



(42) 



The term in square brackets provides a quick way to see the relative importance of the phonoemisive 
process. If we take H^c = lOOOe, then Es^c/fi = 1-8 x lO^^sec^^ We may take D = 0.25 K, and 
i?cx = 5 K as reasonable estimates. For the sound speeds and density, we may take representative values 
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for common organic materials. For paraffin wax, for example, p = 0.91 g/cm'^, cl = 1-94 x 10^ cm/s, and 
we have not been able to find a value for ct- For polyethylene, p = 0.90 g/cm'^, cl — 1.95 x lO'"' cm/s, 
and Ct ~ 0.54 x 10'' cm/s, and for polystyrene, p = 1.06 g/cm'^, cl = 2.35 x 10^ cm/s, and ct = 1.12 x 
10^ cm/s. Hence we take p ~ 1 g/cnv^, and ct = 10^ cm/s. The term in square brackets is then of order 
10~^^, which is disappointingly small. Even if we take A' ~ 100 A as is appropriate if the intermediate 
state is m = ±8 in Feg, the contribution of the phonoemissive process is negligible. Nevertheless, the 
process remains the simplest imaginable candidate for the magnetization relaxation mechanism. To this 
author's mind, it remains an urgent problem to resolve this puzzle. 

Acknowledgements This work was begun with support from the NSF via grant number DMR-0202165. 
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\bm{\epsilon}- 


A 


\bm{\laLnibda} 


Q 


\bm-[\varrho} 




\bm-[\chi} 




\ulll\ \Vd.I ypblXOIlj 


A* 


\ Ulll\ \II1U J 


(T 


\ VmtiJ \ c! t rrm a >■ 
XUlIlL \blgjllld._r 




\Diii\ \psi_r 


t 
s 


\ "hni^ \ "70+" a r- 


ly 


\ "Km -1 \ nn ?• 
\U1i1l \liU.j 


c 
S 


\ "Kn-i \ Trayc! "i crma 

\ UlU I. \ V CLX iD X^UlCLj 




\ "Km -T \ oTTi o era 






? 


\Dml\xij- 










r 


\bin{\itGainina} 




\bm{ \ i t L aunbda} 


i; 


\bm-[\ it Sigma} 




\bm-C\itPsi} 




\bni{\itDelta} 




\bin{\itXi} 


r 


\bm{\itUpsilon} 


n 


\bm{\itOmega} 


@ 


\bin{\itTheta} 


JT 


\bm{\itPi} 




\bm{\itPhi} 







Table 3: Boldface variants of slanted greek letters 



a 


\pmb-[ \upalpha} 


e 


\pmb-[\uptheta} 





\pmb{\upo} 


X 


\pmb-[\uptau} 


P 


\pmb{\upbeta}- 




\pmb-[\upvartheta} 


71 


\pmb-[\uppi} 




\pmb-[ \upups i 1 on} 


Y 


\pmb{ Xupgamma} 


I 


\pmb{\upiota} 


G3 


\pmb{\upvarpi} 




\pmb-[\upphi} 


S 


\pmb{\updelta} 


K 


\pmb{\upkappa} 


P 


\pmb{\uprho} 




\pmb{\upvarphi} 


e 


\pinb{\upepsilon} 


A. 


\pinb{\uplambda} 


Q 


\pmb{\varrho} 


X 


\pinb{\upchi} 


e 


\pmb{\varepsilon} 




\pmb{\upmu} 


a 


\pmb-[\upsigma}- 


V 


\pinb-[\uppsi} 




\pmb-[\upzeta} 


V 


\pmb-[\upnu} 




\pmb{\upvarsigma} 




\pmb{\upomega} 


n 


\pmb{\upeta} 




\pmb{\upxi} 










r 


\bm{\Ganfflia} 


A 


\bm{\Lambda} 


S 


\bin{\Sigma} 




\biii{\Psi} 


A 


\bni{\Delta} 


1— 1 


\biii{\Xi} 


T 


\bni{\Upsilon} 




\biii{\Omega} 





\biii{\Theta} 


n 


\biii{\Pi} 




\biii{\Phi} 







Table 4: Boldface variants of upright greek letters 



2 



